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We study the phase diagram of dense, locally neutral three-flavor quark matter within the frame¬ 
work of the Nambu-Jona-Lasinio model. In the analysis, dynamically generated quark masses are 
taken into account self-consistently. The phase diagram in the plane of temperature and quark che¬ 
mical potential is presented. The results for two qualitatively different regimes, intermediate and 
strong diquark coupling strength, are presented. It is shown that the role of gapless phases diminishes 
with increasing diquark coupling strength. 


I. INTRODUCTION 

Theoretical studies suggest that baryon matter at suf¬ 
ficiently high density and sufficiently low temperature is 
a color superconductor. (For reviews on color supercon¬ 
ductivity see, for example, Ref. [I].) In Nature, the high¬ 
est densities of matter are reached in central regions of 
compact stars. There, the density might be as large as 
lOpo where po ss 0.15 fm -3 is the saturation density. It 
is possible that baryonic matter is deconfined under such 
conditions and, perhaps, it is color-superconducting. 

In compact stars, matter in the bulk is neutral with 
respect to the electric and color charges. Matter should 
also remain in /3 equilibrium. Taking these constraints 
consistently into account may have a strong effect on the 
competition between different phases of deconfined quark 
matter at large baryon densities iiilli . In this pa¬ 
per, we study this competition within the framework of a 
Nambu-Jona-Lasinio (NJL) model. The results are sum¬ 
marized in the phase diagram in the plane of temperature 
(T) and quark chemical potential (p). 

The first attempt to obtain the phase diagram of dense, 
locally neutral three-flavor quark matter as a function of 
the strange quark mass, the quark chemical potential, 
and the temperature was made in Ref. 7J- It was shown 
that, at zero temperature and small values of the strange 
quark mass, the ground state of matter corresponds to 
the color-flavor-locked (CFL) phase @, ! : 9.j. At some crit¬ 
ical value of the strange quark mass, this is replaced by 
the gapless CFL (gCFL) phase .§). In addition, several 
other phases were found at nonzero temperature. For in¬ 
stance, it was shown that there should exist a metallic 
CFL (mCFL) phase, a so-called uSC phase M , as well 


as the standard two-flavor color-superconducting (2SC) 
phase Pi and the gapless 2SC (g2SC) phase [5]. 

In Ref. Q, the effect of the strange quark mass was 
incorporated only approximately through a shift of the 
chemical potential of strange quarks, p s —> p s — m 2 / (2p). 
While such an approach is certainly reliable at small va¬ 
lues of the strange quark mass, it becomes uncontrollable 
with increasing the mass. The phase diagram of Ref. 7| 
was further developed in Refs. iQ where the shift- 
approximation in dealing with the strange quark was not 
employed any more. So far, however, quark masses were 
treated as free parameters, rather than dynamically gen¬ 
erated quantities. In this paper, we study the phase di¬ 
agram of dense, locally neutral three-flavor quark mat¬ 
ter within the NJL model, treating dynamically gener¬ 
ated quark masses self-consistently. Some results within 
this approach at zero temperature were also obtained in 
Refs. UM 

As in Refs. mam, we restrict our analysis to lo¬ 
cally neutral phases only. This automatically excludes, 
for example, mixed 16j] and crystalline [lj] phases. Tak¬ 
ing them into account requires a special treatment which 
is outside the scope of this paper. 

This paper is organized as follows. In Sec. m we in¬ 
troduce the model and, within this model, derive a com¬ 
plete set of gap equations and charge neutrality condi¬ 
tions. The numerical results for the phase diagrams in the 
plane of temperature and quark chemical potential in two 
qualitatively different regimes are presented in Sec. urn 
Finally, our results are summarized in Sec. rm The Ap- 
pendix contains some useful formulas. 
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II. MODEL AND FORMALISM 

In this paper, we use a three-flavor quark model with 
a local NJL-type interaction. The Lagrangian density is 
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given by 

8 
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- K jdet [$ (1 + 75 ) V>] + det [$ (1 - 75 ) i/j\ | , ( 1 ) 

where the quark spinor field ■)/)“ carries color (a = r, g , 6) 
and flavor (a = u, d, s) indices. The matrix of quark 
current masses is given by m = diagy (m u , rrid, m s ). Re¬ 
garding other notations, A a with a = 1, ...,8 are the 
Gell-Mann matrices in flavor space, and Ao = i/2/3 1/. 
The charge conjugate spinors are defined as follows: 
ipc = C'4> T and 1 pc = ^ T C, where ^ = ijj'j 0 is the 
Dirac conjugate spinor and C = * 7 2 7 ° is the charge con¬ 
jugation matrix. 

The model in Eq. 0 should be viewed as an ef¬ 
fective model of strongly interacting matter that cap¬ 
tures at least some key features of QCD dynamics. The 
Lagrangian density contains three different interaction 
terms which are chosen to respect the symmetries of 
QCD. Note that we include the’t Hooft interaction whose 
strength is determined by the coupling constant K. This 
term breaks ?7(1) axial symmetry. 

The term in the second line of Eq. 0 describes a scalar 
diquark interaction in the color antitriplet and flavor an¬ 
titriplet channel. For symmetry reasons there should also 
be a pseudoscalar diquark interaction with the same cou¬ 
pling constant. This term would be important to describe 
Goldstone boson condensation in the CFL phase 0 . In 
this paper, however, we neglect this possibility and there¬ 
fore drop the pseudoscalar diquark term. 

We use the following set of model parameters [Til l: 


m u ,d = 5-5 MeV , (2a) 

m s = 140.7 MeV, (2b) 

G s A 2 = 1.835 , ( 2 c) 

A'A 5 = 12.36 , (2d) 

A = 602.3 MeV . (2e) 


After fixing the masses of the up and down quarks at 
equal values, m U} d = 5.5 MeV, the other four parameters 
are chosen to reproduce the following four observables of 
vacuum QCD 0: m n = 135.0 MeV, m K = 497.7 MeV, 
m v ' = 957.8 MeV, and = 92.4 MeV. This parameter 
set gives = 514.8 MeV 0 . 

In Ref. |13, the diquark coupling Gd was not fixed by 
the fit of the meson spectrum in vacuum. In general, it is 
expected to be of the same order as the quark-antiquark 
coupling Gs- In this paper, we study in detail two pos¬ 
sible cases: the regime of intermediate coupling strength 
with Gd = jGs , and the regime of strong coupling with 
Gd = Gs- The comparison of phase diagrams in these 
two cases will turn out to be very instructive. 


The grand partition function, up to an irrelevant nor¬ 
malization constant, is given by 

Z = e~ nv/T = J VjjVijjSx ( £ +*AW) , (3) 

where f l is the thermodynamic potential density, V is 
the volume of the three-space, and jl is a diagonal ma¬ 
trix of quark chemical potentials. In chemical equilibrium 
(which provides [3 equilibrium as a special case), the non¬ 
trivial components of this matrix are extracted from the 
following relation: 

Mof = (> Q/3 + S a »+[g 3 ( T 3 ) ab + g s (T 8 )J 6°* . 

. , (4) 

Here /1 is the quark chemical potential (by definition, 
g = gs /3 where gs is the baryon chemical potential), 
Hq is the chemical potential of electric charge, while g 3 
and gs are color chemical potentials associated with two 
mutually commuting color charges of the SU(3) C gauge 
group. The explicit form of the electric charge matrix 
is Qf = diag^(|, — A, — A), and the explicit form of the 
color charge matrices is T 3 = diag c (i, — A, 0) and \Z3T S = 
diag c (A, A,-l). 

In order to calculate the mean-field thermodynamic 
potential at temperature T, we first linearize the in¬ 
teraction in the presence of the diquark condensates 
A c ~ (V’c)a* 75 ea/ 3 c£ ofcct/ ; ^ ( n0 sum over c ) an d the quark- 
antiquark condensates a a ~ ( n0 sum over a )■ 

Then, integrating out the quark fields and neglecting the 
fluctuations of composite order parameters, we arrive at 
the following expression for the thermodynamic poten¬ 
tial: 


n — Dl + X l^ c l 2 + 2( - ?s X 

C—1 Q! = l 

T , 5-1 

- 4Ka u a d a s - — X lndet ~jr > (5) 

where we also added the contribution of leptons, f II, 
which will be specified later. 

We should note that we have restricted ourselves to 
held contractions corresponding to the Hartree approx¬ 
imation. In a more complete treatment, among others, 
the’t Hooft interaction term gives also rise to mixed con¬ 
tributions containing both diquark and quark-antiquark 
condensates, i.e., oc cr Q |A a | 2 0. In this study, as 

in Refs. Ij.l2l|. we neglect such terms for simplicity. While 
their presence may change the results quantitatively, one 
does not expect them to modify the qualitative structure 
of the phase diagram. 

In Eq. 0, S^ 1 is the inverse full quark propagator in 
the Nambu-Gorkov representation, 


s - 1 


[G +]" 1 \ 

4>+ [Go]- 1 ) ’ 


( 6 ) 
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with the diagonal elements being the inverse Dirac prop¬ 
agators of quarks and of charge-conjugate quarks, 

[G±]- x = 7 ^±/i7o-M, (7) 

where A' M = (fco,k) denotes the four-momentum of the 
quark. At nonzero temperature, we use the Matsubara 
imaginary time formalism. Therefore, the energy ko is 
replaced with —iui n where ui n = ( 2 n + 1 ) 7 tT are the 
fermionic Matsubara frequencies. Accordingly, the sum 
over K in Eq. © should be interpreted as a sum over 
integer n and an integral over the three-momentum k. 

The constituent quark mass matrix is defined as M = 
diag f (M u , M d , M s ) with 

M a = m a — 4GstJ a + 2A "07307 , ( 8 ) 

where a a are the quark-antiquark condensates, and the 
set of indices (a, (3, 7 ) is a permutation of (u , d , s). 

The off-diagonal components of the propagator © are 
the so-called gap matrices given in terms of three diquark 
condensates. The color-flavor structure of these matrices 
is given by 

M:f=-E^ c ^A c75 , ( 9 ) 

c 

and $ + = 7 °($ _ )l 7 °. Here, as before, a and b refer to 
the color components and a and (3 refer to the flavor com¬ 
ponents. Hence, the gap parameters Ai, A 2 , and A 3 cor¬ 
respond to the down-strange, the up-strange and the up- 
down diquark condensates, respectively. All three of them 
originate from the color-antitriplet, flavor-antitriplet di¬ 
quark pairing channel. For simplicity, the color and flavor 
symmetric condensates are neglected in this study. They 
were shown to be small and not crucial for the qualitative 
understanding of the phase diagram [ 7 ]. 

By making use of the results in Appendix ^ the de¬ 
terminant of the inverse quark propagator can be decom¬ 
posed as follows: 
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( 10 ) 


where €i are eighteen independent positive energy eigen¬ 
values. The Matsubara summation in Eq. © can then 
be done analytically by employing the relation [Hj 


5> 



T 


+ 2 In ^1 + e r 1 ^ . 


( 11 ) 


Then, we arrive at the following mean-field expression for 
the pressure (jp= —H): 


P = 



|e t | + 2 Tln(l + e-Tr) 


+ 4 A'07070,5 - E l Ac ! 2 “ 2Gs E 

c= 1 a=l 

+5s£r" !in ( i+ '' v ) • < i2 > 

l=e,ii e=± ® 


where the contribution of electrons and muons with 
masses m e « 0.511 MeV and ~ 105.66 MeV were 
included. Note that muons may exist in matter in [3 equi¬ 
librium and, therefore, they are included in the model for 
consistency. However, being about 200 times heavier than 
electrons, they do not play a big role in the analysis. 

In this paper, we assume that there are no trapped 
neutrinos in quark matter. This is expected to be a good 
approximation for matter inside a neutron star after the 
short deleptonization period is over. The effect of neu¬ 
trino trapping will be addressed elsewhere {23]. 

The expression for the pressure in Eq. m has a physi¬ 
cal meaning only when the chiral and color superconduct¬ 
ing order parameters, a a and A 0 , satisfy the following set 
of six gap equations: 


dp 

dcr a 

dp 

dK c 


0 , 

0 . 


(13a) 

(13b) 


To enforce the conditions of local charge neutrality in 
dense matter, we also require three other equations to be 
satisfied, 


n Q 

n 3 

n 8 


dp 
dpQ 
dp 
dp 3 
dp 
dp s 


= 0 , 


= 0 


= 0 . 


(14a) 

(14b) 

(14c) 


These fix the values of the three corresponding chemical 
potentials, pq , p 3 and p 8 . After these are fixed, only the 
quark chemical potential p is left as a free parameter. 


III. RESULTS 

In order to obtain the phase diagram, we have to find 
the ground state of matter for each given set of the para¬ 
meters in the model. In the case of locally neutral matter, 
there are two parameters that should be specified: tempe¬ 
rature T and quark chemical potential p. After these are 
fixed, one has to compare the values of the pressure in all 
competing neutral phases of quark matter. The ground 
state corresponds to the phase with the highest pressure. 

Before calculating the pressure, given by Eq. m , one 
has to find the values of the chiral and the color super¬ 
conducting order parameters, a a and A c , as well as the 
values of the three charge chemical potentials, pq , p 3 
and ps■ These are obtained by solving the coupled set of 
six gap equations da together with the three neutrality 
conditions d- By using standard numerical recipes, it 
is not extremely difficult to find a solution to the given 
set of nine nonlinear equations. Complications arise, how¬ 
ever, due to the fact that often the solution is not unique. 

The existence of different solutions to the same set of 
equations, d and d- reflects the physical fact that 
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there could exist several competing neutral phases with 
different physical properties. Among these phases, all but 
one are unstable or metastable. In order to take this into 
account in our study, we look for the solutions of the 
following 8 types: 

1. Normal quark (NQ) phase: Ai = A 2 = A 3 = 0; 

2. 2SC phase: Ai = A 2 = 0, and A 3 7 ^ 0; 

3. 2SCws phase: Ai = A 3 = 0, and A 2 7 ^ 0; 

4. 2SC ds phase: A 2 = A 3 = 0, and Ai 7 ^ 0; 

5. uSC phase: A 2 7 ^ 0, A 3 7 ^ 0, and Ai = 0; 

6 . dSC phase: Ai 7 ^ 0, A 3 7 ^ 0, and A 2 = 0; 

7. sSC phase: A x 7 ^ 0, A 2 7 ^ 0, and A 3 = 0; 

8 . CFL phase: A 3 ± 0, A 2 7 ^ 0, A 3 7 ^ 0. 

Then, we calculate the values of the pressure in all 
nonequivalent phases, and determine the ground state as 
the phase with the highest pressure. After this is done, we 
study additionally the spectrum of low-energy quasipar¬ 
ticles in search for the existence of gapless modes. This 
allows us to refine the specific nature of the ground state. 

In the above definition of the eight phases in terms of 
A c , we have ignored the quark-antiquark condensates a a . 
In fact, in the chiral limit (m a = 0), the quantities cr a 
are good order parameters and we could define additional 
sub-phases characterized by nonvanishing values of one 
or more a a . With the model parameters at hand, how¬ 
ever, chiral symmetry is broken explicitly by the nonzero 
current quark masses, and the values of a a never vanish. 
Hence, in a strict sense it is impossible to define any new 
phases in terms of a a . 

Of course, this does not exclude the possibility of dis¬ 
continuous changes in <j a at some line in the plane of 
temperature and quark chemical potential, thereby con¬ 
stituting a first-order phase transition line. It is generally 
expected that the “would-be” chiral phase transition re¬ 
mains first-order at low temperatures, even for nonzero 
quark masses. Above some critical temperature, however, 
this line could end in a critical endpoint and there is 
only a smooth crossover at higher temperatures. Among 
others, this picture emerges from NJL-model studies, 
both, without |j24 ] and with 0 ] diquark pairing (see also 
Ref. (0). We should therefore expect a similar behavior 
in our analysis. 

Our numerical results for neutral quark matter are 
summarized in Figs. 0 and 0 These are the phase di¬ 
agrams in the plane of temperature and quark chemi¬ 
cal potential, obtained in the mean-field approximation 
in model 0 in the case of an intermediate diquark 
coupling strength, Gd = § Gs, and in the case of a 
strong coupling, Gd = Gs, respectively. The correspond¬ 
ing dynamical quark masses, gap parameters, and three 
charge chemical potentials are displayed in Figs. |o] and 
0 respectively. All quantities are plotted as functions 


of /r for three different fixed values of the temperature: 
T = 0,20,40 MeV in the case of Gd = §Gs (see Fig. 0 
and T = 0,40,60 MeV in the case of Gd = Gs (see 

Fig. 0. 

Let us begin with the results in the case of the diquark 
coupling being Gd = § Gs- In the region of small quark 
chemical potentials and low temperatures, the phase di¬ 
agram is dominated by the normal phase in which the 
approximate chiral symmetry is broken, and in which 
quarks have relatively large constituent masses. This is 
denoted by ySB in Fig. 0 With increasing the tempe¬ 
rature, this phase changes smoothly into the NQ phase 
in which quark masses are relatively small. Because of 
explicit breaking of the chiral symmetry in the model at 
hand, there is no need for a phase transition between the 
two regimes. 

However, as pointed out above, the symmetry argu¬ 
ment does not exclude the possibility of an “accidental” 
(first-order) chiral phase transition. As expected, at lower 
temperatures we find a line of first-order chiral phase 
transitions. It is located within a relatively narrow win¬ 
dow of the quark chemical potentials (336 MeV < ^ < 
368 MeV) which are of the order of the vacuum values 
of the light-quark constituent masses. (For the parame¬ 
ters used in our calculations one obtains M,,. = M,j = 
367.7 MeV and M s = 549.5 MeV in vacuum [H|.) At this 
critical line, the quark chiral condensates, as well as the 
quark constituent masses, change discontinuously. With 
increasing temperature, the size of the discontinuity de¬ 
creases, and the line terminates at the endpoint located 
at (T cr , // cr ) rs (56,336) MeV, see Fig.0 

The location of the critical endpoint is consistent with 
other mean field studies of NJL models with similar 
sets of parameters 0, 0 01. This agreement does not 
need to be exact because, in contrast to the studies in 
Refs. 01110, here we imposed the condition of elec¬ 
tric charge neutrality in quark matter. (Note that the 
color neutrality is satisfied automatically in the normal 
phase.) One may argue, however, that the additional con¬ 
straint of neutrality is unlikely to play a big role in the 
vicinity of the endpoint. 

It is appropriate to mention here that the location of 
the critical endpoint might be affected very much by fluc¬ 
tuations of the composite chiral fields. These are not in¬ 
cluded in the mean-field studies of the NJL model. In 
fact, this is probably the main reason for their inability 
to pin down the location of the critical endpoint consis¬ 
tent, for example, with lattice calculations 0|. (It is fair 
to mention that the current lattice calculations are not 
very reliable at nonzero f .1 either.) Therefore, the predic¬ 
tions of this study, as well as of those in Refs. l2ll 124125| . 
regarding the critical endpoint cannot be considered as 
very reliable. 

When the quark chemical potential exceeds some crit¬ 
ical value and the temperature is not too large, a Cooper 
instability with respect to diquark condensation should 
develop in the system. Without enforcing neutrality, i.e., 
if the chemical potentials of up and down quarks are 
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equal, this happens immediately after the chiral phase 
transition when the density becomes nonzero j2f| . In the 
present model, this is not the case at low temperatures. 

In order to understand this, let us inspect the vari¬ 
ous quantities at T = 0 which are displayed in the up¬ 
per three panels of Fig. 0 At the chiral phase boundary, 
the up and down quark masses become relatively small, 
whereas the strange quark mass experiences only a mod¬ 
erate drop of about 84 MeV induced by the’t Hooft in¬ 
teraction. This is not sufficient to populate any strange 
quark states at the given chemical potential, and the 
system mainly consists of up and down quarks together 
with a small fraction of electrons, see Fig. 0 The elec¬ 
tric charge chemical potential which is needed to main¬ 
tain neutrality in this regime is between about —73 and 
—94 MeV. It turns out that the resulting splitting of the 
up and down quark Fermi momenta is too large for the 
given diquark coupling strength to enable diquark pairing 
and the system stays in the normal phase. 

At h ss 432 MeV, the chemical potential felt by the 
strange quarks, /r —/xq/ 3, reaches the strange quark mass 
and the density of strange quarks becomes nonzero. At 
first, this density is too small to play a sizeable role in 
neutralizing matter, or in enabling strange-nonstrange 
cross-flavor diquark pairing, see Fig. ED The NQ phase be¬ 
comes metastable against the gapless CFL (gCFL) phase 
at /TgCFL ~ 443 MeV. This is the point of a first-order 
phase transition. It is marked by a drop of the strange 
quark mass by about 121 MeV. As a consequence, strange 
quarks become more abundant and pairing gets easier. 
Yet, in the gCFL phase, the strange quark mass is still 
relatively large, and the standard BCS pairing between 
strange and light (i.e., up and down) quarks is not pos¬ 
sible. In contrast to the regular CFL phase, the gCFL 
phase requires a nonzero density of electrons to stay elec¬ 
trically neutral. At T = 0, therefore, one could use the 
value of the electron density as a formal order parameter 
that distinguishes these two phases [6]. 

With increasing the chemical potential further (still at 
T = 0), the strange quark mass decreases and the cross¬ 
flavor Cooper pairing gets stronger. Thus, the gCFL 
phase eventually turns into the regular CFL phase at 
Mcfl ~ 457 MeV. The electron density goes to zero at 
this point, as it should. This is indicated by the vanishing 
value of hq in the CFL phase, see the upper right panel 
in Fig. ED We remind that the CFL phase is neutral be¬ 
cause of having equal number densities of all three quark 
flavors, n u = rid = n s , see Fig. [5] and ED This equality is 
enforced by the pairing mechanism, and this is true even 
when the quark masses are not exactly equal [27^ . 

Let us mention here that the same NJL model at zero 
temperature was studied previously in Ref. 0- Our re¬ 
sults agree qualitatively with those of Ref. @ only when 
the quark chemical potential is larger than the critical 
value for the transition to the CFL phase at 457 MeV. 
The appearance of the gCFL phase for 443 < /r < 
457 MeV was not recognized in Ref. 3|, however. In¬ 
stead, it was suggested that there exists a narrow (about 


12 MeV wide) window of values of the quark chemical po¬ 
tential around /i k, 450 MeV in which the 2SC phase is 
the ground state. By carefully checking the same region, 
we find that the 2SC phase does not appear there. 

This is illustrated in Fig. [ 7 | where the pressure of three 
different solutions is displayed. Had we ignored the gCFL 
solution (thin solid line), the 2SC solution (dashed line) 
would indeed be the most favored one in the interval 
between /i « 445 MeV and /r ss 457 MeV. After including 
the gCFL phase in the analysis, this is no longer the case. 

Now let us turn to the case of nonzero temperature. 
One might suggest that this should be analogous to 
the zero temperature case, except that Cooper pairing 
is somewhat suppressed by thermal effects. In contrast 
to this naive expectation, the thermal distributions of 
quasiparticles together with the local neutrality condi¬ 
tions open qualitatively new possibilities that were ab¬ 
sent at T = 0. As in the case of the two-flavor model 
of Ref. Q, a moderate thermal smearing of mismatched 
Fermi surfaces could increase the probability of creat¬ 
ing zero-momentum Cooper pairs without running into 
a conflict with Pauli blocking. This leads to the appear¬ 
ance of several stable color-superconducting phases that 
could not exist at zero temperatures. 

With increasing the temperature, the first qualita¬ 
tively new feature in the phase diagram appears when 
5 < T < 10 MeV. In this temperature interval, the 
NQ phase is replaced by the uSC phase when the quark 
chemical potential exceeds the critical value of about 
444 MeV. The corresponding transition is a first-order 
phase transition, see Fig. ^ Increasing the chemical po¬ 
tential further by several MeV, the uSC phase is then 
replaced by the gCFL phase, and the gCFL phase later 
turns gradually into the (m)CFL phase. (In this study, 
we do not distinguish between the CFL phase and the 
mCFL phase [t}-) Note that, in the model at hand, the 
transition between the uSC and the gCFL phase is of 
second order in the following two temperature intervals: 
5 < T < 9 MeV and T > 24 MeV. On the other hand, it 
is a first-order transition when 9 < T < 24 MeV. Leav¬ 
ing aside its unusual appearance, this is likely to be an 
“accidental” property in the model for a given set of pa¬ 
rameters. For a larger value of the diquark coupling, in 
particular, such a feature does not appear, see Fig. 0 

The transition from the gCFL to the CFL phase is a 
smooth crossover at all T / 0 00. The reason is that 
the electron density is not a good order parameter that 
could be used to distinguish the gCFL from the CFL 
phase when the temperature is nonzero. This is also con¬ 
firmed by our numerical results for the electric charge 
chemical potential hq i n Fig. ED While at zero tempe¬ 
rature the value of hq vanishes identically in the CFL 
phase, this is not the case at finite temperatures. 

Another new feature in the phase diagram appears 
when the temperature is above about 11 MeV. In this 
case, with increasing the quark chemical potential, the 
Cooper instability happens immediately after the ySB 
phase. The corresponding critical value of the quark che- 
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mical potential is rather low, about 365 MeV. The first 
color superconducting phase is the gapless 2SC (g2SC) 
phase pi- This phase is replaced with the 2SC phase in 
a crossover transition only when /j, > 445 MeV. The 2SC 
is then followed by the gapless uSC (guSC) phase, by the 
uSC phase, by the gCFL phase and, eventually, by the 
CFL phase (see Fig.0. 

In the NJL model at hand, determined by the para¬ 
meters in Eq. <©, we do not find the dSC phase as the 
ground state anywhere in the phase diagram. This is sim¬ 
ilar to the conclusion of Refs. , but differs from that 

of Refs. 30 ■ This should not be surprising because, as 
was noted earlier 0, the appearence of the dSC phase 
is rather sensitive to a specific choice of parameters in 
the NJL model. 

The phase diagram in Fig.0has a very specific ordering 
of quark phases. One might ask if this ordering is robust 
against the modification of the parameters of the model 
at hand. Below we argue that some features are indeed 
quite robust, while others are not. 

It should be clear that the appearance of color¬ 
superconducting phases under the stress of neutrality 
constraints is very sensitive to the strength of diquark 
coupling. In the case of two-flavor quark matter, this was 
demonstrated very clearly in Ref. Q at zero as well as 
at nonzero temperatures. Similar conclusions were also 
reached in the study of three-flavor quark matter at zero 
temperature flSj. 

In the model at hand, it is instructive to study the 
phase diagram in the regime of strong diquark coupling, 
Gd = Gs■ The corresponding results are summarized 
in the diagram in Fig. [2] As we see, the main qualita¬ 
tive difference between the diagrams in Figs. 0 and 0 
occurs at intermediate values of the quark chemical po¬ 
tential. While at Gd = §Gs, there is a large region of 
the g2SC phase sandwiched between the low-temperature 
and high-temperature NQ phases, this is not the case at 
stronger coupling, Gd = Gs- 

The last observation can easily be explained by the 
fact that with increasing diquark coupling strength, the 
condensation energy also increases and therefore Cooper 
pairing is favorable, even if there is a larger mismatch of 
the Fermi surfaces due to charge neutrality constraints. 
Moreover, in the presence of large gaps, the Fermi sur¬ 
faces are smeared over a region of order A. Therefore 
additional thermal smearing is of no further help, and it 
is not surprising that the thermal effects in a model with 
sufficiently strong coupling are qualitatively the same as 
in models without neutrality constraints imposed: ther¬ 
mal fluctuations can only destroy the pairing. In the 
model with a not very strong coupling, on the other hand, 
the interplay of the charge neutrality and thermal fluc¬ 
tuations is more subtle. The normal phase of cold quark 
matter develops a Cooper instability and becomes a color 
superconductor only after a moderate thermal smearing 
of the quark Fermi surfaces is introduced @. 

Other than this, the qualitative features of the phase 
diagrams in Figs. 0 and 0 are similar. Of course, in the 


case of the stronger coupling, the critical lines lie system¬ 
atically at higher values of the temperature and at lower 
values of the quark chemical potential. In this context one 
should note that the first-order phase boundary between 
the two normal regimes “ySB” and “NQ” is insensitive 
to the diquark coupling. Therefore, upon increasing Gd 
it stays at its place until it is eventually “eaten” up by 
the expanding 2SC phase. As a result, there is no longer 
a critical endpoint in Fig. 0 but only a critical point 
where the first-order normal(\SB)-2SC phase boundary 
changes into second order. 


IV. CONCLUSIONS 

In this paper, we studied the T-^jl phase diagram of 
neutral three-flavor quark matter within the NJL model 
of Ref. 31 in which the chiral symmetry is broken explic¬ 
itly by small but nonzero current quark masses. As in the 
previous studies mum , we use the mean-field approx¬ 
imation in the analysis. In contrast to Refs. 1 it flJl Il4| . in 
this paper the constituent quark masses are treated self- 
consistently as dynamically generated quantities. The 
main results are summarized in Figs. 0 and 0 

By comparing our results with those in Ref. 0 (see 
Fig. lib there), we notice several important differences. 
First of all, we observe that a self-consistent treatment 
of quark masses strongly influences the competition be¬ 
tween different quark phases. As was noticed earlier in 
Ref. |3 > there exists a subtle interplay between the two 
main effects. On the one hand, the actual values of the 
quark masses directly influence the competition between 
different normal and color-superconducting phases. On 
the other hand, competing phases themselves determine 
the magnitude of the masses. Very often, this leads to 
first-order phase transitions, in which certain regions in 
the mass-parameter space become inaccessible. 

Some differences to the results in Ref. 0 are related 
to a different choice of model parameters. Most impor¬ 
tantly, the value of the diquark coupling Gd = §Gs is 
considerably weaker than in the NJL model of Ref. 0. 
This can be easily seen by comparing the magnitude of 
the zero-temperature gap at a given value of the quark 
chemical potential, say at /i = 500 MeV, in the two 
models. It is Ag j00 ^ a: 76 MeV in this paper, and it is 

A ( 0 500) « 140 MeV in Ref. Hi. (Note that the strength 
of the diquark pairing in Ref. 11 .‘ill is even weaker, corre¬ 
sponding to Ag 500 -* ss 20 MeV.) It should be noted that 
even the case of the strong coupling, Gd = Gs , which 
corresponds to Aq 500 ' ss 120 MeV, is still slightly weaker 
than that in Ref. 0. In this case, however, the corre¬ 
sponding results differ mostly because the quark masses 
are treated very differently. 

Because of the weaker diquark coupling strength, the 
Cooper instabilities in Fig. 0 happen systematically at 
higher values of the quark chemical potential than in 
Ref. 0. In particular, this is most clearly seen from the 
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critical lines of the transition to the (g)CFL phase. An¬ 
other consequence of the weaker interaction is the pos¬ 
sibility of a thermal enhancement of the (g)2SC Cooper 
pairing at intermediate values of the quark chemical po¬ 
tential. This kind of enhancement was studied in detail 
in Ref. )(|- Making use of the same arguments, one can 
tell immediately how the phase diagram in Fig. ^ should 
change with increasing or decreasing the diquark cou¬ 
pling strength. 

In particular, with increasing (decreasing) the diquark 
coupling strength, the region of the (g)2SC phase at in¬ 
termediate values of the quark chemical potential should 
expand (shrink) along the temperature direction. The re¬ 
gions covered by the other (i.e., uSC and CFL) phases 
should have qualitatively the same shape, but shift to 
lower (higher) values of the quark chemical potential and 
to higher (lower) values of the temperature. In the case 
of strong coupling, in particular, these general arguments 
are confirmed by our numerical calculations. The corre¬ 
sponding phase diagram is shown in Fig. [3 

Several comments are in order regarding the choice of 
the NJL model used here. The model is defined by the 
set of parameters in Eq. o which were fitted to repro¬ 
duce several important QCD properties in vacuum IS]. 
(Note that the same model was also used in Ref. Q.) 
It is expected, therefore, that this is a reasonable effec¬ 
tive model of QCD that captures the main features of 
both chiral and color-superconducting pairing dynamics. 
Also, a relatively small value of the cutoff parameter in 
the model, see Eq. l!2ell . should not necessarily be viewed 
as a bad feature of the model. In fact, this might simply 
mimic a natural property of the full theory in which the 
coupling strength of relevant interactions is quenched at 
large momenta. 

In this relation, note that the approach of Ref. ^3 re¬ 
garding the cutoff parameter in the NJL model is very 
different. It is said there that a large value of this pa¬ 
rameter is beneficial in order to extract results which are 
insensitive to a specific choice of the cutoff. However, we 
do not find any physical argument that would support 
this requirement. Instead, we insist on having an effective 
model that describes reasonably well the QCD properties 
at zero quark chemical potential. We do not pretend, of 
course, that a naive extrapolation of the model to large 
densities can be rigorously justified. In absence of a better 
alternative, however, this seems to be the only sensible 
choice. 

The results of this paper might be relevant for under¬ 
standing the physics of (hybrid) neutron stars with quark 
cores, in which the deleptonization is completed. In order 
to obtain a phase diagram that could be applied to pro¬ 
toneutron stars, one has to generalize the analysis to take 
into account neutrino trapping. This work is in progress 
now 

In the end, it might be appropriate to mention that, 
despite the progress in our understanding of the phase 
diagram of neutral dense quark matter, there still exists 
a fundamental problem here. The reason is that some 


regions of the phase diagrams in Figs.n anc 10 corres P on d 
to phases that are known to be unstable ( 23 . Of course, 
it is of prime importance to resolve this issue. 

Note added. While writing our paper, we learned that a 
partially overlapping study is being done by D. Blaschke, 
S. Fredriksson, H. Grigorian, A.M. Ozta§, and F. Sandin 

Ill- 


Acknowledgments 


This work was supported in part by the Virtual 
Institute of the Helmholtz Association under grant 
No. VH-VI-041 and by Gesellschaft fiir Schwerionen- 
forscliung (GSI) and by Bundesministerium fiir Bildung 
und Forschung (BMBF). S. R. thanks for using the Cen¬ 
ter for Scientific Computing (CSC) of the Johann Wolf¬ 
gang Goethe-Universitat Frankfurt am Main. 


APPENDIX A: EVALUATION OF THE 
THERMODYNAMIC POTENTIAL 

We use the following ordering of the quark field com¬ 
ponents: 

■ (Ai) 

In this basis, the matrices of quark current and con¬ 
stituent masses read 


m = diag (m„, m d , m s , m u , m d , m s , m u , m d , m s ) (A2) 
and 

M = diag (M u , M d , M s , M u , M d , M s , M u , M d , M s ) , 

(A3) 

respectively, with M a given by Eq. ©■ Moreover, the ma¬ 
trix of quark chemical potentials takes the general form 


£ = diag(<,^,^,/x»,^,/xf,^,/x^,^) . (A4) 


Finally, the explicit color-flavor structure of the gap ma¬ 
trices 4 >± [see Eq. ©] is given by 


/ 


$ = -75 


V 


0 0 0 0 
0 0 0 -A: 
0 0 0 0 
0 -a 3 0 0 

A 3 0 0 0 
0 0 0 0 
0 0 -a 2 0 
0 0 0 0 
A 2 0 0 0 


A 3 0 0 0 
0 0 0 0 
0 0 -a 2 0 
0 0 0 0 
0 0 0 0 
0 0 0 —Ai 

0 0 0 0 
0-Ai0 0 
Ai 0 0 0 


a 2 \ 
0 
0 
0 

Ai 

0 

0 

0 

0 ) 


(AS) 

and $ + = — (4> - )^. By making use of the symmetries in 
the model, we choose all three parameters A c to be real. 




In order to calculate the last term in the thermody¬ 
namical potential in Eq. it is useful to rewrite the 
determinant of the inverse full propagator as det ,!? -1 = 
det (7 o7o 5' _1 ) = det( 7 0 )det (yo-S 1-1 ) = det ( 70 S' -1 ). 
One should also note that the matrices 7 o[Cj ] _1 and 
7 o < f >± can be expressed in terms of the spin projectors, 


7o[G 0 ± ]- 1 = £ 

S 

7 o $ ± = ±^2 


k 0 ± /1 — M —sk 

—sk kg ± fi + M 

0 $ 

0 


Vs , (A 6 ) 
(A7) 


where the two projectors are defined as 


1 


V s = -(1 + s cr ■ k) , for s = ±. (A 8 ) 


Here k = k /k, k = |k|, and $ represents only the 
color-flavor part of the gap matrix in Eq. El. i.e., 
= ± 75 $. By making use of the definition in Eq. 
as well as Eqs. ifAlill and E3, we obtain the following 
representation: 


70 S-^Ys^Vs, 


(A9) 


where S s 

_ 

ito 

— M s , and 



/- 

M + 

M sk 

0 

M s = 


sk 

—fi — M 

-<l 


0 

-<% 

fi + M 


l 

$ 

0 

sk 

(with s = 

±) 

is real and symmetric. S 


$ 

0 

sk 

fi — M J 


(A10) 


plicit energy dependence in M s , their eigenvalues e,; de¬ 
termine the quasiparticle dispersion relations, k 0 = ti{k). 
By using the matrix relation lndet(A) = Trln(A) as well 
as the properties of projectors V s , we derive 

In det (yoS 1-1 ) = In ^detS 1 ^ 1 • det^Z 1 ^ . (All) 


It turns out that the two determinants appearing on the 
right hand side of this equation are equal, i.e., det SZ = 
det^jl 1 . From the physics viewpoint, this identity re¬ 
flects the twofold spin degeneracy of the spectrum of 
quark quasiparticles. The formal proof of this degener¬ 
acy is straightforward after noticing that the following 
matrix relation is satisfied: 


szl = ns^K- 1 , (A 12 ) 


where 


Another observation, which turns out to be helpful in 
the calculation, is that the determinant det>§Z 1 (fco) is 
an even function of fco, he., detSg (—to) = det S^Bto)- 
This is a formal consequence of the following matrix re¬ 
lation: 

Sr 1 (-* 0 ) = -BS 7 1 (to) tf -1 . (ah) 

where the explicit form of the unitary matrix B is 


/ 0 0 0 i\ 

0 0 -i 0 I 

0 i 0 0 

\-i 0 0 0/ 


(A15) 


It satisfies det ,8 = 1. The invariance of the determinant 
det (to) with respect to the change of the energy sign, 
k 0 —> —ko, is directly related to the use of the Nambu- 
Gorkov basis for quark fields. In this basis, for each quasi¬ 
particle excitation with a positive energy kg = e(k), there 
exists a corresponding excitation with a negative energy 
to = —e(k). Therefore, the result for the determinant 
should read 

18 

det (S'” 1 ) = (hi - e 2 i) 2 . (A16) 

i=1 

In order to simplify the numerical calculation of the 
eigenvalues of the matrix M+, defined in Eq. (fTTfl , we 
first write it in a block-diagonal form. The total dimen¬ 
sion of this matrix is 36 x 36. With a proper ordering 
of its rows and columns, it decomposes into 6 diagonal 
blocks of dimension 4x4 and one diagonal block of di¬ 
mension 12 x 12. The explicit form of these blocks reads 


M { 1 ] = 


( -Vd + M d k 

k - M d 


0 

V —a 3 


a 3 

0 


0 

a 3 

Mu + M u 


—a 3 
0 
k 


k fi- M u 

(Al7a) 


Mf = 


(fi - M d k 0 A 3 

k Md + M d A 3 0 

0 A 3 -fi - M u k 

V -A 3 0 k -fi + M u , 

(A17b) 


Mf = 


( —fi +M S k 

k -fi - M s 


0 

^2 


-A 2 \ 

0 


0 


A 2 Mu + M u 


\ -a 2 


0 


k Mu - M u ) 


(Al7c) 


/l 

0 

0 

-1 

0 

0 

0 

(A13) 


/mS 

-M 

k 

0 

—a 2 \ 

0 

0 

-1 

0 

m (4) - 


k 

Ml + M s 

a 2 

0 

Vo 

0 

0 

1) 


JVI + — 


0 

a 2 

—Mu - M u 

k 






l 

-A 2 

0 

k 

_ Mu + M u J 


is a unitary matrix with unit determinant, det 1Z = 1. (A17d) 
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and 


M^ = 


l-fi 9 s + M s k 0 

k -/zf - M s Ai 
0 Ai fi d + M d 

\ -A x 0 k ■■ b 


-A, \ 

0 


(A17e) 


Mf = 


k 

0 

V - A i 


Ms 


0 

Ai 


“Ai \ 
0 


Ai ~n b d - M d k 
0 k —/z^ "d M d J 

(Al7f) 


J 


/—MU - M u k 

k -/£ + M, 


M 7) = 


0 

0 

0 

0 

0 

—a 3 

0 

0 

0 

—a 2 


0 

0 

0 

0 

A 3 

0 

0 

0 

A 2 

0 


0 
0 

Mu - M U 

k /if 
0 

a 3 
0 
0 

0 

A 2 
0 
0 


■ 


-A 3 

0 

0 

0 

—a 2 
0 
0 
0 


0 

0 

0 

—a 3 

~i4l - M d 
k 
0 
0 

0 

0 

0 

-Ai 


Out of 36 eigenvalues from all seven blocks, there are 
18 positive and 18 negative eigenvalues. Out of total 18 
positive eigenvalues, 9 of them correspond to quark type 
quasiparticles and the other 9 correspond to antiquark 
type quasiparticles. In our calculation, we extract all 36 
eigenvalues numerically and then use them in the calcu¬ 
lation of the pressure, see Eqs. m through 

Here, it might be interesting to note that the eigenval¬ 
ues of the 4x4 matrices can be calculated analytically 
in the limit when two quark masses appearing in each of 
them are equal. For example, when M d = M u , the four 
eigenvalues of matrix + 1 are given by 





± s/Ml + k^j +A2 


Hd- V 9 u 

2 

(A18) 


0 0 -A 3 

0 a 3 0 

A 3 0 0 

0 0 0 

k 0 0 

+ M d 0 0 

0 H 9 d -M d k 
0 k n 9 d + M d 

0 0 -Ai 

0 Ai 0 

Ai 0 0 

0 0 0 


0 0 

0 0 

0 A 2 

—a 2 0 

0 0 

0 0 

0 Ai 

-A 3 0 

~H b - M s k 

k -ft + M s 

0 0 [i\ 

0 0 


0 -A 2 \ 

a 2 0 

0 0 

0 0 

0 -Ai 

Ai 0 

0 0 

0 0 

0 0 

0 0 

-M a k 

k ^ b + M s J 

(A17g) 


( 2 ) 

while the eigenvalues of Aid differ only by the sign in 
front of the second term, 


Af , = ± ^(ri+^ ±v ^ J|TF ) 2 + al + ifAA. 

(A19) 

When the value of SM = M d — M u is nonzero but 
small, the c orrections to the above eigenvalues are 
±M U SM / (2^/M% + k 2 ) with the plus sign in the case 
of antiparticle modes, and the minus sign in the case of 

particle modes. The eigenvalues of A4 < | 3 ' and A4+* in the 

(5) 

limit M s = M u , as well as the eigenvalues of Add and 
Mf in the limit M s = M d , are similar. 
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FIG. 1: The phase diagram of neutral quark matter in the regime of intermediate diquark coupling strength, Gd = § Gs- 
First-order phase boundaries are indicated by bold solid lines, whereas the thin solid lines mark second-order phase boundaries 
between two phases which differ by one or more nonzero diquark condensates. The dashed lines indicate the (dis-)appearance 
of gapless modes in different phases, and they do not correspond to phase transitions. 



FIG. 2: The phase diagram of neutral quark matter in the regime of strong diquark coupling, Gd — Gs- The meaning of the 
various line types is the same as in Fig. Q 
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FIG. 4: The dependence of the quark masses, of the gap parameters and of the charge chemical potentials on the quark chemical 
potential at a fixed temperature, T = 0 MeV (three upper panels), T = 40 MeV (three middle panels), and T = 60 MeV (three 
lower panels). The diquark coupling strength is Gd = Gs • 
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FIG. 5: The dependence of the number densities of quarks and electrons on the quark chemical potential at T = 0 MeV for 
diquark coupling strength Gd = |Gs- Note that the densities of all three quark flavors coincide above /x = 457 MeV. The 
density of muons vanishes for all values of /x. 



FIG. 6: The dependence of the number densities of quarks, electrons and muons on the quark chemical potential at T = 0 MeV 
for diquark coupling strength Gd = Gs■ Note that the densities of all three quark flavors coincide above /x = 414 MeV. 
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FIG. 7: Pressure divided by /i 4 for different neutral solutions of the gap equations at T = 0 as functions of the quark chemical 
potential /r: regular CFL (bold solid line), gapless CFL (thin solid line), 2SC (dashed line), normal (dotted line). The diquark 
coupling strength is Gd = §Gs- 
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